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> ■ ABSTRACT 

l> 
(N 

q \ A key parameter to the description of all star formation processes is the 

■<^|- \ density structure of the gas. In this letter, we make use of probability distri- 

bution functions (PDFs) of Herschel column density maps of Orion B, Aquila, 
and Polaris, obtained with the Herschel Gould Belt survey (HGBS). We aim to 
understand which physical processes influence the PDF shape, and with which 
signatures. The PDFs of Orion B (Aquila) show a lognormal distribution for low 
column densities until Ay ~ 3 (6), and a power- law tail for high column densities, 
consistent with a p oc r~ 2 profile for the equivalent spherical density distribution. 
The PDF of Orion B is broadened by external compression due to the nearby 
OB stellar aggregates. The PDF of a quiescent subregion of the non-star-forming 
Polaris cloud is nearly lognormal, indicating that supersonic turbulence governs 
the density distribution. But we also observe a deviation from the lognormal 
shape at Ay>l for a subregion in Polaris that includes a prominent filament. 
We conclude that (i) the point where the PDF deviates from the lognormal form 
does not trace a universal Ay -threshold for star formation, (ii) statistical density 
fluctuations, intermittency and magnetic fields can cause excess from the lognor- 
mal PDF at an early cloud formation stage, (iii) core formation and/or global 
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collapse of filaments and a non-isothermal gas distribution lead to a power-law 
tail, and (iv) external compression broadens the column density PDF, consistent 
with numerical simulations. 



1. Introduction 

The star formation process represents a dramatic transformation of a molecular cloud 
in time and space where the main governing elements are turbulence, gravity, and magnetic 
fields. The spatial structure of clouds, now impre ssively revealed by H ersche l imaging ob- 
servations in the far-infrared (e.g., Andre et al. ( 120101 ). Motte et al. (120 lOh . Molinari et 



al. (120101 )). is ver y inh omogeneous and dominated by filaments (Arzoumanian et al. (120111) . 
Schneider et al. ( 120121 )). It is only with Herschel space observatory (Pilbratt et al. (120101 )) 
observations that diffuse to dense gas is n ow tr aced at high angular resoluti on (ty pically 
18"). Combining PACS (Poglitsch et al. (bold )) and SPIRE (Griffin et al. (bold )) data 
provides column density maps that are superior to tho se ob tained from extinction usin g 
near-IR data (Lombardi et al. ( 120061 ). Kainulainen et al. (120091 ). Froebrich & Rowles ( 120101 ). 
Schneider et al. (120 llh ) that have angular resolutions of ~2' and suffer from saturation at 
visual extinctions Ay above ^25. 

Here, we aim to disentangle the relative contributions of turbulence, gravity, and exter- 
nal compression that influence the density structure of a molecular cloud. A useful analysis 
technique is to use probability distribution functions (PDFs) of the column density, which 
characterizes the fract ion of gas with a column density N in the range [N, N+AN] (e.g., 
Federrath et al. (120101 )). Extinction maps (see above) have shown that molecular clouds can 
have a lognormal PDF for low column densities, and either a power-law tail or more complex 
shapes for higher column densities. Is otherm al, hydrodynamic simulations including turbu- 
lence and gravity (e.g., Klessen et al. (J2000h ) have shown that gravitational collapse i nduce s 
a power-law tail in th e PD F at high densities. More recent studies (Kritsuk et al. ( 120111 ). 
Federrath & Klessen (120 131 ) and references therein) have investigated which parameters in- 
fluence the shape of the PDF. Following these studies, fitting the slope of the high-density 
tail of the PDF allows us to determine the exponent a of an equivalent spherical density 
distribution p{r) = p (r/r )~ a . 

In this study we make use of Hersciel-derived column density PDFs of the Orion B 
molecular cloud, a template region for studies of low- to high-mass star formation (Lada 



1 Herschel is an ESA space observatory with science instruments provided by European- led Principal 
Investigator consortia and with important participation from NASA. 
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et al. (1l99lh ) . Orion B is amongst the nearest (distance ~400 pc; Gibb (120081 )) giant 
molecular cloud (GMC) complexes, with a mass of around 10 5 M , and hosts several OB- 
clusters (NGC2023/24, NGC2068/71). Orion B is located within the Ha-shell 'Barnard's 
Loop' and diverse OB stellar aggregates impact the cloud from the west with radiation and 
stellar winds. To understand better what governs the density structure and its link to star 
formation, we compare the Orion B PDFs to those obtained with the HGBS for a quiescent 
cloud (Polaris) and a star- forming region (Aquila). 



2. Observations 



Orion B, Aquila, and Polaris were observed with the PACS and SPIRE in strum ents 
onboard Herschel as part of the Herschel Gould Belt survey (HGBS, Andre et al. ( 120101 )) in 
parallel mode with a scanning speed of 60"/sec and two orthogonal coverages. The Orion B 
data w ere ob tained on 2010 Sep 29 and 2011 M ar 13. For details on Pola ris se e Men'shchikov 
et al. (boioh . Miville-Deschen es et a l. (J201oh . Ward-Thompson et al. fcoioh . and Konyves 
et al. ( )2010l ). Bontemps et al. ( I2010l ) for Aquila. The angular resolutions at 160 /tm (PACS), 
250 /im, 350 /im, and 500 /mi (all SPIRE), are ~12", ~18", ~25", and ~36", respectively. 
The SPIRE data were reduced with HIPE version 7.1956, including a destriper-module with 
a polynomial baseline of zeroth order. Both scan directions were then combined using the 
'naive-mapper', i.e., a simple averaging algorithm. The PACS data were reduced using 
HIPE 6.0.2106. In addition to the standard data reduction steps, non-linearity correction 
was applied on the 160 /im signal, which affects only the bright (> 10 Jy /pixel) regime. The 
levell data were then combined into a map with Scanamorphos vlO (Roussel (120121 )). 



Column density and dust temperature maps were determined fro m a m odified blackbody 
fit to the wavelengths 160 to 500 /mi (see, e.g., Konyves et al. (120 10h ). We recovered 
the Herschel zer o-flux levels of the Orion B field for each wavelength with Planck data 
(Bernard et al. (J2010l )). For the region covered by both PACS and SPIRE simultaneously, 
we fixed the specific dust opacity per unit mass (dust+gas) appro ximated by the power law 
k u = 0.1 (y J 1000G H Zf cm 2 /g and (3=2 (cf. Hildebrand (119831 )). took a mean molecular 
weight per hydrogen molecule of 2.8, and left the dust temperature and column density as 
free parameters. As an im provement to this procedure, we applied the technique described 
in Palmeirim et al. ( (120131 )) that uses the flux information of the 500 /mi map but with the 
help of a multi-resolution decomposition, the angular resolution of the final maps is higher, 
i.e., that of the 250 /im data at 18". To test the robustness of the derived high- resolution 
map of Orion B, we constructed ratio maps between the 18"-resolution column density map 
- smoothed to the common-resolution of the 500 /mi data (36")- and the originally 36"- 
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resolution maps. This ratio map has a mean value of 1.0 and a standard deviation of 0.03. 
The two column density maps agree within 15%. In ad dition , we investigated the effect of 



increasing opacity for high column densities (Roy et al. (( 120131 )) on the PDF and found that 
its dispersion decreases (~10-20%) and can provoke a steeper slope of the power-law tail for 
high densities. 

Because the density structure of molecular clouds depends on how energy is injected 
into a cloud (spiral density waves, expanding supernovae shells, H n-regions, or gravitational 
contraction), we determine the hydrodynamic Mach-number Ai that characterizes to first or- 
der the influence of turbulence. Stronger isothermal, non-magnetized supersonic turbulence 
leads to a higher Mach number and thus stronger local density enha ncements. In contrast, 



magnetic fields smooth out density variations (Molina et al. (120121 )). M. can be derived 
from observations of the full width at half maximum (FWHM in [km s -1 ]) of a molecular 
line, and the sound speed c s =0.188^T kin /10K with the kinetic temperature T fcin . 

M = (V3FWHM)/(c s V81n2) (1) 

If the LTE assumption is valid, Tu n ~T e:E , with T ex = 5.53[ln(5.53/T m b) + l)]" 1 for the 
optically thick 12 CO 1— >0 line. If gas and dust are well mixed, the temperature should also 
correspond to the dust temperature derived from Herschel. We emphasize, however, that the 
determination of the Mach-number remains rather uncertain (error ~30-40%) and mainly 
gives a tendency. 



3. The column density structure of Orion B 

3.1. Column density maps 

The column density map of Orion B (Fig. [1]) is dominated by the two active star- 
forming clumps NGC2023/24 and NGC2068/71 (Buckle et al. 2010) with very high local 
column densities N(H 2 ) up to a few 10 22 cm -2 and high dust temperatures of up to 35 K 
(Fig. [2]) due to the H n-regions. These two dense ridges are outlined by a column density 
level of ~3-4xl0 21 cm -2 and stand out in an extended cloud with a typical column density 
of l-2xl0 21 cm -2 . In contrast, the northern-eastern part of Orion B is colder and less 
active. A sharp cut-off in column density at the western border of NGC2023/2024 and the 
Horesehead nebula is seen in column d ensity cuts at constant declination (Fig. [I]) which was 



known from CO data (Wilson et al. (120051 )). From west to east, we first observe a strong 



increase of column density on a few pc scale for NGC2023/2024, and a weaker but clearly 
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visible increase for the southern region. East of the peaks, the column density decreases to 
a level of 1 to 2xl0 21 cm" 2 which is higher than the v alues a t the western border. Such a 



profile was also seen in the Pipe nebula (Peretto et al. ( 120121 )) where the authors proposed 
a large-scale compression by the winds of the Sco OB2 association likely caused this sharp 
edge. A similar process may be at work for Orion B because the western part of the cloud 
is exposed to various OB aggregates (OBlb-d). 



3.2. Probability distribution functions of column density 

The distribution of number of pixels vs. column density for Orion B, Aquila, and Polaris 
are displayed in Fig. HJ We will use the term probability distribution function (PDF) and 
the notation p(r/) though the pure pixel distribution is not strictly a PDF which is defined 
for a lognormal distribution as 

p(rj)drj = (2 n a 2 )" ' 5 exp[-(77 - /i) 2 /(2a 2 )] drj (2) 

with rj = ln(N/{N)) and <J n as the dimensionless dispersion of the logarithmic field, and [i 
the mean. The normalization allows a direct comparison between clouds of different column 
density, and a v is a measure for the density variation in a turbulent medium. We determine 
a v with a fit to the assumed lognormal low-density part of the PDF and the slope with the 
index s from a power-law fit with p{rj) = po{v/Vo) s to the high-density part. The results of 
the PDF fit are given together with the Mach- number determination in Table 1. 

The Orion B and Aquila PDFs show a well-defined lognormal part for low column 
densities and a clear power-law tail at higher column densities, starting at an extinction^ 
Ay around 3 for Orion B and 6 for Aquila. To first order, the PDFs only differ in their 
width (0^=0.45 for Orion B and 0.3 for A quila) . The PDFs of two subregions (Fig. [3]) in 



the Polaris cirrus cloud (Falgarone et al. (119981 )) are more narrow with a=0.22 and 0.27 



respectively. The PDF of the quiescent region is almost perfectly lognormal, however, above 
Ay~1.5 we observe a slight excess which ist most likely a resolution effect but may possibly 
result from a physical process (see Sec. 3.3). A clear deviation from the lognormal shape is 
found for the 'saxophone' filament. Interestingly, the excess for Ay >1 has not the form of a 
power-law tail. Note that all PDF features (shape, width etc.) do not depend on the angular 
resolution (18" vs. 36") or pixel number (our statistic here is high because the images are on 
a 6" grid). To quantify better the deviation of the PDF from lognormal, we determined the 



2 



For better comparison to the literature values, we use the visual extinction value derived f rom th e column 



density adopting the conversion formula N(H 2 )/Av=0.94xl0 21 cm 2 mag 1 (Bohlin et al. (119781 )) 
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Table 1: Temperature regime, Mach-number, and PDF fit results for Orion B, Aquila, and 
Polaris 
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(1) Observed excit ation temperature range f rom 1 2 CO 1— >0 data (Ori on B: Buckle et al. ( 2010h . Aquila: 
Zeilik et al. |l978h : Polaris: Bcnsch ct al. (J2003I ). Shimoikura et al. (J2012h ). No large-scale CO data is 
availabe for Polaris and Aquila). 

(2) Observed dust temperature range from Herschel data. 

(3) and (4) average temperature. 

(5) line width from 12 CO l->0. 

(6) sonic Mach-number from average temperature and CO-line width. 

(7) dispersion of the PDF. 

(8) exponent of the spherical density profile. 
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higher momentqj skewness S and kurtosis /C. The skewness, describing the asymmetry of the 
distribution, is positive for Orion B, Aquila, and Polaris-saxophone (5=1.17, 1.24, and 0.49, 
error typically 0.05), implying an excess at higher column densities, and is 5=0.19 (/C=3.0) 
for Polaris-quiet, confirming the nearly lognormal form of its PDF. The Orion B and Aquila 
PDFs have much higher values for the kurtosis (X=6.8 and 7.9) that arise from pronounced 
wings. The Polaris-saxophone region has also a (high) value of /C=4.4, indicating an excess 
at high column densities. 

From theory, a purely lognormal distribution is only expected if the cloud structure is 
shaped by supersonic, isothermal turbulence while deviations in the form of a po wer law 



for high column densities are p redicted for self-gravitating clouds (Klessen et al. ( 120001 ) 



Ballesteros-Paredes et al. fj201l[ )). The concept of isothermality does not fully apply to all 
clouds, as can be seen in the temperature PDFs in Fig. [2J While Polaris can be considered 
as a nearly isothermal gas phase, Orion B and Aqu ila sh ow a more complex temperature 



distribution over a larger range. Froebrich & Rowles (j2010l ) argued that the Ay-value where 



the transition of the PDF must take place is alw ays a round 6, and defined this value as a 



threshold for star formation. Kainulainen et al. (120111 ) found values between Ay=2-5 and 
proposed a scenario in which this Ay-range marks a transition between dense clumps and 
cores and a more diffuse interclump medium. The PDFs of Orion B and Aquila shown in 
Fig. H] clearly show that the transition from lognormal to power law is not universal but varies 
between clouds (Ay ~3 and 6, respectively). This behaviour suggests that the Ay -transition 
value neither represents a universal threshold in star formation nor a phase transition (unless 
the density of the clumps and the interclump medium strongly varies from cloud to cloud) . 
A similar result was obtained in the study of Schneider et al. (2013), presenting a large 
sample of PDFs from low- to high-mass star-forming clouds. Deviations of the PDF from 
lognormal are more likely a function of cloud parameters, in particular the virial parameter, 
the dominant forcing model, and the Mach number as shown in models (Federrath & Klessen 
2013) and observations (Schneider et al. 2013). 

From the slope of the power-law tail for Orion B and Aquila, the exponent of an equiv- 
alent spherical density distribution p(r) oc r~ a is determined to be a=1.99 (1.77) for Orion 
B (Aquila), conforming with results typically obtained for individual collapsing cores. The 
high-density tail, however, cannot be explai ned by the core population alone because it does 



not provide sufficient mass (Konyves et al. (j2010l )). Moreover, it is probably also caused by 



3 <S = -^ j 00 ^ dr]p{ri)[r) - (tj)] 3 and /C = -^ J^ drip(r])[r] - (r/)] 4 (see, e.g., Federrath et al. 2010) 

4 Compressive modes on large to small scales are generated by galactic spiral shocks, expanding super- 
nova shells and Hll regions, gravitational contraction, and outflows. Solenoidal forcing arises from galactic 
rotation and magneto-rotational instabilities. 



global grav itation al collapse of large r spatial areas like fil ament s and ridges (see, e.g., Schnei- 
der et al. (120 10[ ). Hill et al. (120 111 ). Palmeirim et al. (120 131 ) for observational examples). 
It was shown that non-isothermal flows can also cause power-law tails (Passot &Vazquez- 
Semadeni 1998) so that this process may influence the shape of the PDF as well. Indeed, 
in regions with significant temperature variations, a can reach values lar ger th an free-fall 
(a=2.4 for the high-mass star-forming cloud NGC6334, Russeil et al. (120131 )). possibly 
pointing towards a scenario in which heating/cooling processes become important. 



3.3. Comparison to models 

We compare our observational PDFs with those obtained from hydrodynamic simula- 
tions (Federrath & Klessen 2013), including gravity, magnetic fields, and different turbulent 
states (.M=2-50, star formation efficiencies from to 20%, and different forcing modes). In 
addition, we determine the ' b-para meter' (a 2 = f 2 a 2 = In (1 + b 2 A4 2 ), e.g., Federrath et al. 
2010, Burkhart & Lazarian (120121 )). characterizing the link between density and velocity in 
a cloud. We find: 

(1) The dispersion a v of the PDF for Aquila is 0.3 and for the two Polaris subregions 0.22 
and 0.27, while cr r; =0.45 for Orion B. At the same time, the Mach number for all regions is 
typically 6 to 8 (in view of the uncertainty of Ai, they are bascially the same), while only 
the Polaris-quiet subregion has a significant lower value of 3. Numerical models indicate 
that a larger width is caused by a higher Mac h num ber and/or comp ressiv e forcing instead 
of solenoidal forcing (see also Federrath et al. ( 120101 ). Tremblin et al. (120 121 ) ). Since Orion B 
and Aquila have similar values of Mach number, we conclude that the Orion B cloud is likely 
exposed to compressive modes - as seen also in the sharp cutoff of column density (Fig. [T]) 
- caused by the stellar winds from diverse OB aggregate^. Aqu ila has been proposed to 
be located at an encounter of several superbubbles (Frisch (1l998h ). but the impact on its 
density structure - more or less important than close-by OB-stars - can not b e inferred. 
Both clouds are exposed to relatively high magnetic fields (Crutcher et al. ( 119991 ). Sugitani 
et al. ( 1201 ll )). so that the more narrow PDF of Aquila is presumably not caused by magnetic 
fields alone. 

(2) Polaris-quiet has a narrow (<r=0.22), lognormal PDF, the gas is nearly isothermal (see 
temperature PDF in Fig. [2]), and has a low (~3) Mach- number. Only isothermal turbulence 
simulations without self-gravity reproduce this shape of the PDF. We computed the forcing- 



parameter b = 1/M. x (exp ((/ a Ti 



\0.5 



using an average of 2.5 between solenoidal and 



5 Note that the Pipe cloud as a clear example of compression (Sec. 3.1) also shows a broad PDF with 
r^=0.60 in extinction maps (Schneider et al. 2013.) 
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compressive forcing (Federrath et al. 2010) for the factor / = cr s /a v (estimation of the 3D 
density fluctuation a s out of the 2D column density fluctuation a n ). The resulting value of 
6=0.2 is lower than what was found by compari ng; wi th the purely solenoidal driven isother- 



mal MHD simulations of Burkhart & Lazarian (120 12l )). Our data point for Polaris-quiet fits 
on their model (Fig. 3) with b=l/3. In any case, these results show that the Polaris-quiet 
PDF is consistent with the view that the cloud's density distribution is mainly governed by 
solenoidal forcing. 

3) Power-law tails in the high-density PDF regime form under the presence of self-gravity, 
but ca n also be provoked by purely non-isothermal turbulence (Passot & Vazquez-Semadeni 



(119981 )). For Orion B and Aquila, gravity most likely dominates because a large number of 
pre- and protostellar dense cores and supercritical filaments are present (Andre et al. 2010, 
Konyves et al., in prep.). The gas is not isothermal (FigHJ), but the temperature does not 
vary by several orders of magnitude either. The excess in the PDF for the Polaris-saxophone 
region is more difficult to interpret. In this gravitational an d the rmally subcritical fil amen t, 



only one candidate pre-stellar core (Ward-Thompson et al. (J2010J), Shimoikura et al. (120121 )) 
was found. The gas can be condidered as isothermal, so that here magnetic fields may play 
a role (the strength is not known), leading to a narrow PDF (e.g., Molina et al. (2012)), 
or statistical density fluctuations and intermittency due to locally compressive turbulence. 
These effects may also explain the slight excess in the PDF of Polaris-quiet. 

SPIRE has been developed by a consortium of institutes led by Cardiff Univ. (UK) and 
including: Univ. Lethbridge (Canada); NAOC (China); CEA, LAM (France); IFSI, Univ. 
Padua (Italy); IAC (Spain); Stockholm Observatory (Sweden); Imperial College London, 
RAL, UCL-MSSL, UKATC, Univ. Sussex (UK); and Caltech, JPL, NHSC, Univ. Colorado 
(USA). This development has been supported by national funding agencies: CSA (Canada); 
NAOC (China); CEA, CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB (Swe- 
den); STFC, UKSA (UK); and NASA (USA). We thank J.-Ph. Bernard for providing the 
flux values from Planck. C. Federrath acknowledges the Australian Research Council for a 
Discovery Projects Fellowship (grant No. DP110102191). Part of this work was supported 
by the ANR-11-BS56-010 project "STARFICH". 
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Orion B H 3 — column density 



i? 

6T)0 00 



cut in [arcmin] 
100 50 







7 10 s1 



6 10? 1 



10" 



5 55 00 



5 50 00' 



5 45 m 00 a 
RA(3000) 



5 40 00 



5"35 m 00 B N (h 2 ) 

[cm~ ; 



Fig. 1. — H 2 -column density map at 18" angular resolution of Orion B obtained from Herschel 
data. Known H n regions are labeled. The panel inside the image shows cuts (color-coded in 
blue, green, and red) in H2-column density of the NGC2023/24 region at constant declination. 
These cuts are indicated in the image. 
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Orion B dust temperature 
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Fig. 2. — Dust temperature map at 18" angular resolution of Orion B obtained from Herschel 
data. The panel inside the image shows the temperature PDFs of Orion B, Aquila, and 
Polaris (the whole region, not seperated by subregions). 
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Fig. 3. — Hersehel SPIRE map at 250 yum of the Polaris and Aquila fields. The two regions 
for which we show the corresponding column density PDFs for Polaris are indicated by 
dashed white lines. These are the 'saxophone', a region containing a filament with high 
column density but lacking star formation and a subregion further north-east. 
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Fig. 4. — Probability distribution functions of column density for Orion B and Aquila (top) 
and Polaris (bottom) at an angular res olutio n of 18" (a similar version of the Aquila PDF 
was previously shown in Andre et al. ( 120111 )). The left y-axis is the PDF as a normalized 
probability p(rj) while the right y-axis indicates the number of pixels per logarithmic bin. 
Note that due to the large number of pixels (6" grid), the error-bars calculated using Poisson- 
statistics are very small. The PDFs do not change using a lower sampled grid (i.e. at the 
resolution of 18" or 36"). The lower x-axis gives the column density N(H2) in units of 10 21 
cm -2 (corresponding approximately to the Ay in magnitudes using the Bohlin-factor) . The 
upper x-axis is the dimensionless parameter 77 = \n(N/(N)). The green curve indicates the 
fitted lognormal PDF and the red line the power law fit to the high-density tail. The width 
of the fitted lognormal PDF (cr,,), the powerx-law slope index s, and the exponent a of 
the equivalent spherical density profile (p oc r~ a ) are given in each panel. Note that the 
variation of opacity with density causes a systematic error on the PDF which is larger than 
the statistical error. We estimate the error on a and a to be typically 10-20%. 



